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We show that, in a system with defects, two-particle states may experience destructive quantum 
interference, or antiresonance. It prevents an excitation localized on a defect from decaying even 
where the decay is allowed by energy conservation. The system studied is a qubit chain or an 
equivalent spin chain with an anisotropic (XXZ) exchange coupling in a magnetic field. The chain 
has a defect with an excess on-site energy. It corresponds to a qubit with the level spacing different 
from other qubits. We show that, because of the interaction between excitations, a single defect 
may lead to multiple localized states. The energy spectra and localization lengths are found for two- 
excitation states. The localization of excitations facilitates the operation of a quantum computer. 
Analytical results for strongly anisotropic coupling are confirmed by numerical studies. 
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I. INTRODUCTION 

One of the most important potential applications of 
quantum computers (QC's) is studies of quantum many- 
body effects. It is particularly interesting to find new 
many-body effects in condensed-matter systems that 
could be easily simulated on a QC. In the present paper 
we discuss one such effect: antiresonance, or destructive 
quantum interference between two-particle excitations in 
a system with defects. We also study interaction-induced 
two-particle localization on a defect and discuss implica- 
tions of the results for quantum computing. 

The basic elements of a QC, qubits, are two-state sys- 
tems. They are naturally modeled by spin-1/2 particles. 
In many suggested realizations of QC's, the qubit-qubit 
interaction is "on" all the time [lj. In terms of spins, it 
corresponds to exchange interaction. The dynamics of 
such QC's and spin systems in solids have many impor- 
tant similar aspects that can be studied together. 

In most proposed QC's the energy difference between 
the qubit states is large compared to the qubit-qubit in- 
teraction. This corresponds to a system of spins in a 
strong external magnetic field. However, in contrast to 
ideal spin systems, level spacings of different qubits can 
be different. A major advantageous feature of QC's is 
that the qubit energies can be often individually con- 
trolled 0, 0, 3 This corresponds to controllable disor- 
der of a spin system, and it allows one to use QC's for 
studying a fundamentally important problem of how the 
spin-spin interaction affects spin dynamics in the pres- 
ence of disorder. 

Several models of QC's where the interqubit inter- 
action is permanently "on" are currently studied. In 
these models the effective spin-spin interaction is usu- 
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ally strongly anisotropic. It varies from the essentially 
Ising coupling <r*crLin nuclear magnetic resonance and 
some other systems [MBSSIyl {n,m enumerate qubits, 
and z is the direction of the magnetic field) to the 
XF-type (i.e., a^a^ + cr^m) or the XXZ-type (i.e., 
AnmofjOf^ -(- a^a^ + ^(T^) coupling in some Josephson- 
j unction based systems @,|3j]. 

The Ising coupling describes the system in the case 
where the transition frequencies of different qubits are 
strongly different. Then cr^a^ is the only part of the 
interaction that slowly oscillates in time, in the Heisen- 
berg representation, and therefore is not averaged out. If 
the qubit frequencies are close to each other, the terms 
a n a m + <J n <J m become smooth functions of time as well. 
They lead to resonant excitation hopping between qubits. 
In a multiqubit system with close frequencies, both Ising 
and A"F interactions are present in the general case, but 
their strengths may be different In this sense the 

XXZ coupling is most general, at least for qubits with 
high transition frequencies. 

The interqubit interaction often rapidly falls off with 
the distance and can be approximated by nearest neigh- 
bor coupling. Many important results on anisotropic spin 
systems with such coupling have been obtained using the 
Bethe ansatz. Initiall y th e emphasis was placed on sys- 
tems without defects [lj] or with defects on the edge of 
a spin chain |l3| ■ More recently these studies have been 
extended to systems with defects that are described by 
integrable Hamiltonians [1J|. However, the problem of 
a spin chain with several coupled excitations and with 
defects of a general type has not been solved. 

In this paper we investigate interacting excitations in 
an anisotropic spin system with defects. We show that 
the excitation localized on a defect does not decay even 
where the decay is allowed by energy conservation. We 
also find that, in addition to a single-particle excitation, 
a defect leads to the onset of two types of localized two- 
particle excitations. 

The analysis is done for a system with the XXZ cou- 
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pling. The coupling anisotropy is assumed to be strong, 
as in the case of a QC based on electrons on helium, 
for example 0|. The ground state of the system corre- 
sponds to all spins pointing in the same direction (down- 
wards, for concreteness) . A single-particle excitation cor- 
responds to one qubit being excited, or one spin being 
flipped. If the qubit energies are tuned in resonance with 
each other, a QC behaves as an ideal spin system with 
no disorder. A single-particle excitation is then magnon- 
type, it freely propagates through the system. 

In the opposite case where the qubit energies are tuned 
far away from each other (as for diagonal disorder in 
tight-binding models), all single-particle excitations are 
localized. If the excitation density is high, the interac- 
tion between them may affect their localization, leading 
to quantum chaos, cf. Refs. llCM. Il5l ll6tll71 Understand- 
ing the interplay between interaction and disorder is a 
prerequisite for building a QC. We will consider the case 
where the excitation density is low, yet the interaction 
is important. In particular, excitations may form bound 
pairs (but the pair density is small). 

One of the important questions is whether the interac- 
tion leads to derealization of excitations. More specifi- 
cally, consider an excitation, which is localized on a de- 
fect in the absence of other excitations. We now create an 
extended magnon-type excitation (a propagating wave), 
that can be scattered off the localized one. The problem 
is whether this will cause the excitation to move away 
from the defect. We show below that, due to unexpected 
destructive quantum interference, the scattering does not 
lead to derealization. 



A. Model and preview 

We consider a one-dimensional array of qubits which 
models a spin-1/2 chain. For nearest-neighbor coupling, 
the Hamiltonian is 



Jii (J n <J n+l 



(1) 
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Here, a\ are the Pauli matrices and U = 1. The parame- 
ter J characterizes the strength of the exchange coupling, 
and A determines the coupling anisotropy. We assume 
that |A| » 1; for a QC based on electrons on helium J A| 
lies between 20 and 8, for typical parameter values Q. 

We will consider effects due to a single defect. Respec- 
tively, all on-site spin-flip energies are assumed to be 
the same except for the site n = hq where the defect is 
located, that is, 



:■(«) 



(2) 



In order to formulate the problem of interaction- 
induced decay of localized excitations, we preview in 
Fig. a part of the results on the energy spectrum of 
the system. In the absence of the defect, the energies 



of single-spin excitations (magnons) lie within the band 
El ± J, where e\ = e — J A (the energy is counted off 
from the ground-state energy) . The defect has a spin- flip 
energy that differs by g (a qubit with a transition fre- 
quency different from that of other qubits). It leads to 
a localized single-spin excitation with no threshold in g 1 
for an infinite chain. The energy of the localized state is 
shown by a dashed line on the left panel of Fig. ^ 
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FIG. 1: Left panel: the one-excitation energy spectrum in an 
infinite spin chain with a defect. The energies of extended 
states (magnons) form a band of width J centered at E\ . The 
dashed line indicates the energy of the excitation localized on 
the defect. Right panel: the two-excitation energy spectrum. 
The band 2(e\ ± J) is formed by uncoupled magnons. The 
band centered at 2ei + g is formed by localized-delocalized 
pairs (LDP's) in which one excitation is localized on the de- 
fect and the other is in an extended state. The narrow band 
centered at 2ei + J A is formed by propagating bound pairs 
(BP's) of excitations. The dashed lines show the energies of 
the states where both excitations are bound to the defect. 

We now discuss excitations that correspond to two 
flipped spins. A defect-free XXZ system has a two- 
magnon band of independently propagating noninteract- 
ing magnons. However, the anisotropy of the exchange 
coupling leads also to the onset of bound pairs (BP's) 
of excitations. The BP band is much narrower than the 
two-magnon band and is separated from it by a compar- 
atively large energy difference JA, see the right panel 
of Fig. n In the presence of a defect, there are two- 
excitation states with one excitation localized on the de- 
fect and the other being in an extended state. We call 
them localized-delocalized pairs (LDP's). An interplay 
between disorder and interaction may lead to new types 
of states where both excitations are localized near the 
defect. Their energies are shown in the right panel of 
Fig-fflby dashed lines. 

A localized one-spin excitation cannot decay by emit- 
ting a magnon, by energy conservation. But it might 
experience an induced decay when a magnon is inelas- 
tically scattered off the excited defect into an extended 
many-spin state. Magnon-induced decay is allowed by 
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energy conservation when the total energy of the local- 
ized one-spin excitation and the magnon coincides with 
the energy of another two-particle state. In the XXZ 
model the total number of excitations (flipped spins) is 
conserved, and therefore decay is only possible into ex- 
tended states of two bound magnons. In other words, it 
may only happen when the LDP band overlaps with the 
BP band in Fig. □ 

Decay into BP states may occur directly or via the 
two-excitation state located next to the defect. The am- 
plitudes of the corresponding transitions turn out to be 
nearly equal and opposite in sign. As a result of this 
quantum interference, even though the band of bound 
magnons is narrow and has high density of states, the 
LDP to BP scattering does not happen, i.e., the excita- 
tion on the defect is not delocalized. The BP to LDP 
scattering does not happen either, i.e., a localized exci- 
tation is not created as a result of BP decay. 

In Sec. II below and in the Appendix we briefly ana- 
lyze localization of one excitation in a finite chain with a 
defect, for different boundary conditions. In Sec. Ill wc 
discuss the two-excitation states localized near a defect. 
In Sec. IV we consider the resonant situation where the 
energy band of extended bound pair states is within the 
band of energies of the flipped defect spin plus a magnon, 
i.e., where the BP band overlaps with the LDP band in 
Fig. ^ We find that the localized excitation remains on 
the defect site in this case. Analytical results for a chain 
with strong anisotropy |A| are compared with numerical 
calculations. Section V contains concluding remarks. 



their localization are interesting for quantum computing. 
Coherent excitation transitions allow one to have a QC 
geometry where remote "working" qubits are connected 
by chains of "auxiliary" qubits, which form "transmission 
lines" |^. Localization, on the other hand, allows one to 
perform single-qubit operations on targeted qubit. 

A QC makes it possible to model spin chains with dif- 
ferent boundary conditions. The simplest models are an 
open spin chain with free boundaries, which is mimicked 
by a finite-length array of qubits (for electrons on helium, 
it can be implemented using an array of equally spaced 
electrodes, cf. Ref. \19t) or a periodic chain, which can be 
mimicked by a ring of qubits. 

An open A-spin chain is described by the Hamiltonian 
Q), where the first sum runs over n = 1, . . . , N and the 
second sum runs over n = 1, . . . , N — 1 (the edge spins 
have neighbors only inside the chain). In what follows, we 
count energy off the ground-state energy E = —(Ne + 
g)/2 + (A - I)JA/4, i.e., we replace in Eq. Q H -> 
H-E . 

The eigenfunctions of H in the case of one excitation 
can be written as 0] 

N 

ip 1 = J2a(n)4>(n), (4) 

n=l 

where 4>(n) corresponds to the spin on site n being up and 
all other spins being down. The Schrodinger equation for 
a(n) has the form 



II. ONE EXCITATION: LOCALIZED AND 
EXTENDED STATES 

In order to set the scene for the analysis of the two- 
excitation case, in this section and in the Appendix 
we briefly discuss the well-known case of one excitation 
(flipped spin) in an XXZ spin chain with a defect 0] 
and the role of boundary conditions. The Hamiltonian 
of the chain with the defect on site no has the form 
©. We assume that the excitation energy e largely ex- 
ceeds both the coupling constant | J| and the energy ex- 
cess on the defect site \g\. In this case the ground-state of 
the system corresponds to all spins being parallel, with 
(<7^) = — 1 irrespective of the signs of J, g, and A. 

Without a defect, one-spin excitations are magnons. 
They freely propagate throughout the chain. The term 
in the Hamiltonian Q responsible for one-excitation hop- 

ping is iZ (t) = J2 n H n \ with 

HrP = jJj2i=x,y a n a n+l = + (3) 

Ar) _ r.Wlf _ _+ - 
En — [In \ — cr„ +1 (7 n . 

(r) (I) 

The operators t n and t n cause excitation shifts n — > 
n + 1 and n + 1 — > n, respectively. 

A defect leads to magnon scattering and to the onset 
of localized states. Both propagation of excitations and 



/ . JA r JA r \ . . 

I £l + gOn,n + —°ri,l + —0 n ,N I a(n) 

+ -[a(n- 1) + a(n+ 1)] = E x a{n), (5) 

where E\ — e — J A is the energy of a flipped spin in an 
ideal infinite chain in the absence of excitation hopping 
and E\ is the one-excitation energy eigenvalue. For an 
open TV-spin chain we set a(0) = a(N + 1) = in Eq. JSJ. 

The Hamiltonian of a closed A-spin chain has the form 
where both sums over n go from I to A and the 
site A + I coincides with the site 1. Here, the defect 
location no can be chosen arbitrarily. The wave func- 
tion can be sought in the form |0J. The Schrodinger 
equation then has the form (JSJ), except that there are no 
terms proportional to <$ n ,i) 3n,N from the end points of 
the chain. It has to be solved with the boundary condi- 
tion a(n + A) = a(n). 

For an open chain, the solution of the Schrodinger 
equation JSJ can be sought in the form of plane waves 
propagating between the chain boundaries and the de- 
fect, 

o(n) = C\ r e l9n + Cl r e- Wn , |n - no| > 1. (6) 

The subscripts I and r refer to the coefficients for the 
waves to the left (n < no) and to the right (n > no) from 
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the defect. The interrelations between these coefficients 
and the coefficient a(n ) follow from the boundary con- 
ditions and from matching the solutions at no- They are 
given by Eqs. (|A1|) and [|A2(1 . 

For a closed chain, on the other hand, the solution can 
be sought in the form 

a(n) = Ce l9n + C 'e~ ien ; n <n<N, (7) 
a[n) = C e mn+N ^ +C'e- l0{n+N \ l<n<n . 

The energy E\ as a function of 9 can be obtained by 
substituting Eqs. © and into Eq. Q. Both for the 
open and closed chains it has the form 

Ei = E X + J cos 9. (8) 

The eigenfunctions ©, Q) with real 9 correspond to 
sinusoidal waves (extended states). From JHJ, their en- 
ergies lie within the band £\ ± J. In contrast, local- 
ized states have complex 9, and their energies lie outside 
this band. The corresponding solutions for both types of 
chains are discussed in the Appendix. In a sufficiently 
long chain, there is always one localized one-excitation 
state on a defect. Its energy is given by Eq. (|A8|) . In 
the case of an open chain with the coupling anisotropy 
parameter |A| > 1, there are also localized states on the 
chain boundaries. Their energy is given by Eq. (|A6|> . 

III. TWO EXCITATIONS: UNBOUND, BOUND, 
AND LOCALIZED STATES 

A spin chain with a defect displays rich behavior in 
the presence of two excitations. It is determined by the 



interplay between disorder and inter-excitation coupling. 
Solutions of the two-excitation problem have been ob- 
tained in the case of a disorder potential of several spe- 
cial types, where the system is integrable 14]. Here we 
will study the presumably nonintegrable but physically 
interesting problem where the on-site energy of the defect 
differs from that of the host sites. 

The system is described by Hamiltonian JQ). In or- 
der to concentrate on the effects of disorder rather than 
boundaries, we will consider a closed chain of length 
N. We will also assume that the anisotropy is strong, 
|A| » 1. 

The wave function of a chain with two excitations is 
given by a linear superposition 



tp2 = /J a{n, m)(j){n, m), (9) 

n<m 



where 4>{n, to) is the state where spins on the sites n 
and to are pointing upward, whereas all other spins are 
pointing downward. In a periodic chain, the sites with 
numbers that differ by N are identical, therefore we have 
a(n, m) = a(m, n + N). 

From Eq. the Schrodinger equation for the coeffi- 
cients a(n, to) is 



(2si + gS n 
1 



JAS, 



m,n+l 



a(n, to) 



+ - J[a(n - 1, to) + [a(n + 1, to) + a(n, to — !)](! - S m , n+1 ) + a(n, to + 1)] 



E2a{n, to). 



(10) 



r 



Here, E^ is the energy of a two-excitation state \e\ — 
e — J A is the on-site one-excitation energy, cf. Eq. ©]. 
As before, we assume that the defect is located on site 
n . 



A. An ideal chain 

In the absence of a defect the system is integrable. The 
solution of Eq. ljTt)|> can be found using the Bethe ansatz 

a(n,m) = Ce t{9in+e27 ^ + C e i (^ n +^ m l (11) 



The energy of the state with given 9i , 92 is obtained by 
substituting Eq. i|ll[l into Eq. I|1(J|I written for to > n + 1. 
This gives 



E 2 = 2s i + J(cos6»i + cos 9 2 ). 



(12) 



By requiring that the ansatz l|ll|) apply also for to = 
n + 1, we obtain an interrelation between the coefficients 
C and C, 



C_ 

C 7 



1 - 2Ae l<>1 + e l ^ 1+e ^ 
1 - 2Ae i02 + e^ 8l+s ^ ' 



(13) 



With account taken of normalization, Eqs. (|llfl and i|13|) 
fully determine the wave function. The states with real 
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6*i 2 form a two-magnon band with width 4|J|, as seen 
from the dispersion relation l|12l) . The magnons are not 
bound to each other and propagate independently. 

For |A| > 1, Eq. (fbl)) also has a solution C — 0, which 
gives a complex phase 6*2 = 0* = |$ — it with k > 0. 
This solution corresponds to the wave function a(n, m) oc 
exp[i9(n + m)/2 — n(m — n)]; we have m > n. From 
Eqs. (JTTJl-lO we obtain 



J 



e- K = A" 1 cos(0/2), E BP = Eg}, + ^ costf, 



Eg'p = 2e x + J A - 



2A' 



(14) 



Eq. I|14|) describes a bound pair of excitations. Such 
a pair can freely propagate along the chain. The wave 
function is maximal when the excitations are on neigh- 
boring sites. The size of the BP, i.e. the typical distance 
between the excitations, is determined by the reciprocal 
decrement « , and ultimately by the anisotropy param- 
eter A. For large |A|, the excitations in a BP are nearly 
completely bound to nearest sites. Then the coefficients 
a(n,m) oc 6 n +i imi to the lowest order in |A| _1 . 

The distance between the centers of the BP band and 
the two-magnon band Epp — 2e\ is given approximately 
by the BP binding energy JA. The width of the BP 
band | J/A| is parametrically smaller than the width of 
the two-magnon band 4|J|, see Fig.^ 

For nearest-neighbor coupling and for large |A|, trans- 
port of bound pairs can be visualized as occurring via an 
intermediate step. First, one of the excitations in the pair 
makes a virtual transition to the neighboring empty site, 
and as a result the parallel spins in the pair are separated 
by one site. The corresponding state differs in energy by 
J A from the bound-pair state. At the next step the sec- 
ond spin can move next to the first, and then the whole 
pair moves by one site. From perturbation theory, the 
bandwidth should be J 2 /J A = J/ A, which agrees with 
Eq. (T2J. 

The above arguments can be made quantitative by in- 
troducing an effective Hamiltonian = Yln=i of 
BP's. It is obtained from the Hamiltonian H^> © in the 
second order of perturbation theory in the one-excitation 
hopping constant J, 



H (t) 



J 
64A 

J 
64A 



u(r) Al) Ar) i 

I'ti- l L n— 1 ' L n+l L n+l\ 

Y u n L n-1 ' L n L ji+l\- 



(15) 



The operators tn'^ of excitation hopping to the right 
or left are given by Eq. 1(5} . The first pair of terms in 

H$p [Eq. I|15l) ] describes virtual transitions in which a 
BP dissociates and then recombines on the same site. 
This leads to a shift of the on-site energy level of the BP 
by J/2 A. The second pair of terms describes the motion 
of a BP as a whole to the left or to the right. 



The action of the Hamiltonian on the wave func- 
tion a(n, n + 1) is given by 



H {t} a(n,n+ 1) 
J 



(16) 



J 



2A 



a(n, n + 1) + ^^I a ( n + 1 ; n + 2) + a(n — 1, n)]. 



The Schrodinger equation for a BP is given by the sum 
of the diagonal part (the first term) of Eq. l|10fl with m = 
n + 1 and the right-hand side of Eq. (|16(l . In this approxi- 
mation a BP eigenfunction is <z(n, m) = 5 n +\ t m exp(i9n), 
and the dispersion law is of the form (|14|) . 



B. Localized states in a chain with a defect 

We now consider excitations in the presence of a defect. 
In this section we assume that the defect excess energy g 
is such that 



| 5 |»|J|, |JA- 5 |»|J|. 



(17) 



The first inequality guarantees that the localization 
length of an excitation on the defect is small [its inverse 
Im 9 d « In2|#/J| > 1 cf. Eq. 

The second condition in Eq. (|17J) can be understood 
by noticing that, in a chain with a defect, there is a two- 
excitation state where one excitation is localized on the 
defect whereas the other is in an extended magnon-type 
state. These excitations are not bound together. The en- 
ergy of such an unbound localized-delocalized pair (LDP) 
should differ from the energy of unbound pair of magnons 
by « g, and from the energy of a bound pair of magnons 
()14J) by w J A — g. In this section we consider the case 
where both these energy differences largely exceed the 
magnon bandwidth J (the case where | JA — g\ < | J| will 
be discussed in the following section). 



1. Unbound localized-delocalized pairs (LDP's) 

In the neglect of excitation hopping, the energy of an 
excitation pair where one excitation is far from the de- 
fect (| n — no | ^ 1) and the other is localized on the 
defect is 2ei + g. [This can be seen from Eq. (|10|) for the 
coefficients a(n, m) in which off-diagonal terms are dis- 
regarded.] At the same time, if one excitation is on the 
defect and the other is on the neighboring site no±l, this 
energy becomes 2ei + g + J A. The energy difference JA 
largely exceeds the characteristic bandwidth J. There- 
fore, if the excitation was initially far from the defect, it 
will be reflected before it reaches the site tiq ± 1. 

The above arguments suggest to seek the solution for 
the wave function of an LDP in a periodic chain in the 
form 



a(n ,m) = Ce l8m + C*'e^ m , 



(18) 



with the boundary condition a(riQ, no + 1) = a(no, tiq + 
N — 1) = 0. This boundary condition and the form of 
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the solution are similar to what was used in the problem 
of one excitation in an open chain. 

From Eq. {It)]), the energy of an LDP is 



Eldp ~ 2ei + g + J cost 



(19) 



The wave number 9 takes on N — 3 values 7rfc/ (N — 2), 
with k — 1, . . . , N — 3. As expected, the bandwidth of 
the LDP is 2\J\, as in the case of one-excitation band in 
an ideal chain. 

We have compared Eq. (|19fl with numerical results ob- 
tained by direct solution of the eigenvalue problem (|10|l . 
For N = 100, A = 20, g/J = 10 we obtained excellent 
agreement once we took into account that the energy 
levels lfH?|) are additionally shifted by J 2 /2g. This shift 
can be readily obtained from Eq. (|10fl as the second-order 
correction (in J) to the energy Ed of the excitation lo- 
calized on the defect. It follows also from Eq. (|A8|) for 

\g\ » \J\- 

We note that the result is trivially generalized to the 
case of a finite but small density of magnons. The 
wave function a(no, mi, 7712, • • ■ , mu) of M uncoupled 
magnons and an excitation on the site no is given by 
a sum of the appropriately weighted permutations of 
exp(±i#imi ± i#2?7i2 ■ ■ ■±*^m'TIm) over the site numbers 
rrii (these numbers can be arranged so that mi < . . . < 
ink < no < mk+i < . . . < TOm)> with real 9i. The weight- 
ing factors for large |A| are found from the boundary 
conditions and the condition that a(no,mi,m,2, . . .) =0 
whenever any two numbers rrii,mi-\-i differ by 1. When 
the ratio of the number of excitations M to the chain 
length N is small, the energy is just a sum of the energies 
of uncoupled magnons and the localized excitation, i.e., 

it is [M + l)ei +g + J J2iLi cos For sma11 M / N : mat- 
tering of a magnon by the excitation on a defect occurs 
as if there were no other magnons, i.e., the probability of 
a three-particle collision is negligibly small. 



2. Bound pairs localized on the defect: the doublet 

For large |A|,|p/J|, a bound pair of neighboring ex- 
citations should be strongly localized when one of the 
excitations is on the defect site. Indeed, if we disregard 
intersite excitation hopping, the energy of a BP sitting 

on the defect is E^ — 2si + g + J A. It differs signifi- 
cantly from the energy of freely propagating BP's (|14f> . 
causing localization. 

The major effect of the excitation hopping is that the 
pair can make resonant transitions between the sites 
(no, no + 1) and (no — l,no)- Such transitions lead to 
splitting of the energy level of the pair into a doublet. 
To second order in J the energies of the resulting sym- 
metric and antisymmetric states are 

(Q), J(2JA + 9)± ^ 



E 



D 



E 



D 



4A( J A + g) 4( J A + g) 



The energy splitting between the states is small if J A and 
g have same sign. If on the other hand, | J A + g\ < \J\, 



the theory has to be modified. Here, the bound pairs 
with one excitation localized on the defect are resonantly 
mixed with extended unbound two-magnon states. We 
do not consider this case in the present paper. 

We note that, in terms of quantum computing, the on- 
set of a doublet suggests a simple way of creating entan- 
gled states. Indeed, by applying a ir pulse at frequency 
£\ + g one can selectively excite the qubit no. If then 
a 7r pulse is applied at the frequency — (si + g), it 
will selectively excite a Bell state [ 1 1 ) + |10)] /y/2, where 
\ij) describes the state where the qubits on sites n — 1 
and no + 1 are in the states \i) and \j), respectively. The 
excitations on sites no ± 1 can then be separated without 
breaking their entanglement using two-qubit gate opera- 
tions. 



3. Localized states split off the bound-pair band 

The energy difference E^ — Ebp ~ g between BP 
states on the defect site l|20[l and extended BP states 
largely exceeds the bandwidth \J/A\ of the extended 
states. Therefore it is a good approximation to assume 
that the wave functions of extended BP states are equal 
to zero on the defect. In other words, such BP's are 
reflected before they reach the defect. In this sense, the 
defect acts as a boundary for them. One may expect that 
there is a surface-type state associated with this bound- 
ary. 

The emergence of the surface-type state is facilitated 
by the defect-induced change of the on-site energy of a 
BP located next to the defect on the sites (no + 1, no + 2) 
[or (no — 2, no — 1)]. This change arises because virtual 
dissociation of a BP with one excitation hopping onto a 
defect site gives a different energy denominator compared 
to the case where a virtual transition is made onto a 
regular site. It is described by an extra term SEbp in the 
expression (1 1511 for the diagonal part of the Hamiltonian 

n + 1 [or n = n — 2] , 



Hn^ with n 



5E BP = 



Jg 



4A(JA-g) 



(21) 



Using the transformed Hamiltonian (|16f) . one can an- 
alyze BP states in a way similar to the analysis of one- 
excitation states in an open chain, see the Appendix. The 
Schrodinger equation for BP states away from the defect 
is given by the first term in Eq. I jlOj l and by Eqs. i|16|) 
and (|21|l . The BP wave functions can be sought in the 
form 



a(n,n+ 1) = Ce Wn + C'e 



(22) 



Then the BP energy as a function of 9 is given by Ebp 



(o) 

BP 



(J/2 A) cos 6», cf. Eq. O. 



E 

The values of 9 can be found from the boundary condi- 
tion that the BP wave function is equal to zero on the de- 
fect, i.e., a(no, no + 1) = a(no — 1, no) = a(no+N— 1, no + 
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N) = 0. From the Schrodinger equations for a(n, n + 1) 
with n = n + 1 and n = N + n — 2, we obtain an 
equation for 9 of the form 

2 



/(0) = /(-0), /(0) = 



5E BP - ^ 



J8(N-3) 



(23) 



Equation has 2(N - 1) solutions for exp(i0). The 
solutions exp(£0) = ±1 are spurious, in the general case. 
The roots 9 and —9 describe one and the same wave func- 
tion. Therefore there are N — 2 physically distinct roots 
9, as expected for an TV-spin chain with two excluded BP 
states located at (n ,n ± 1). 

Depending on the ratio q = J/(4A5Ebp) = (JA — 
g) / g, the roots are either all real or there is one or two 
pairs of complex roots with opposite signs. Real roots 
correspond to extended states, whereas complex roots 
correspond to the states that decay away from the de- 
fect. The onset of complex roots can be analyzed in the 
same way as described in the Appendix for one excita- 
tion. There is much similarity, formally and physically, 
between the onset of localized BP states next to the de- 
fect and the onset of surface states at the edge of an open 
chain. 

We rewrite Eq. i|23|) in the form 



tan[0(JV - 1)] = 



2 sin 9 (cos 9 — q) 



(24) 



{cos 9 -q) 2 -sin 2 9 

For \q\ > 1 all roots of Eq. l(2*4"f are real. At \q\ = 1 there 
occurs a bifurcation where two real roots with opposite 
signs merge (at 9 — 0, for q = 1, or at 9 = n, for q = —1). 
They become complex for | q | < 1. For \q\ = 1— 2(N— 1) _1 
two other real roots coalesce at 9 = or it and another 
pair of complex roots emerges. 

As the length of the chain increases, the difference be- 
tween the pairs of complex roots decreases. In the limit 
N — > oo the roots merge pairwise. The imaginary part 
of one of the roots is 



Im 9 BP — \n\SE BP /(J/4A)\. 



(25) 

The second root has opposite sign. 

One can show from the Schrodinger equation for the 
BP's that the solution 125fl describes a "surface- type" BP 
state localized next to the defect on sites (no + 1, no + 2). 

(s) 

The solution with — 9 BP describes the surface state on 
(no — 2, no — 1). The amplitudes of these states expo- 
nentially decay away from the defect. For J A > g > 
JA/2 > or JA < g < JA/2 < we have Re 0^j> = 0. 
For g > J A > or g < J A < 0, we have Re 9 pp = n, 
and decay of the localized state is accompanied by oscilla- 
tions. The complex roots in a finite chain can be pictured 
as describing those same states on the opposite sites of 
the defect. But now the states are "tunnel" split because 
of the overlap of their tails inside the chain, which leads 
to the onset of two slightly different localization lengths. 
The energy of the localized state in a long chain is 



E 



M 

BP 



E 



(0) 
BP 



SE 



BP 



W4A) ; 
SEpp 



It lies outside the band E^ p ± ( J/2A) JS) of the ex- 
tended BP states. The distance to the band edge strongly 
depends on the interrelation between the defect excess 
energy g and the BP binding energy JA. It is of the 
order of the BP bandwidth | J/A|, except for the range 
where the difference between JA and g becomes small. In 
this range the energy of the surface-type state sharply in- 
creases in the absolute value. This is illustrated in Fig. [21 
The localization length |Im 9^ p \~ 1 is large when the state 
energy is close to the band edge and shrinks down with 
decreasing | J A — g\, i.e., with increasing \SEbp\- 




1 1.5 

g/(JA) 



FIG. 2: The distance between the energy level of the BP local- 
ized next to the defect and the center of the band of extended 
BP states, e B p = {E^l - £^,)/(J/4A), vs the scaled defect 
excess energy g/JA. The localized state exists in an infinite 
chain for g/ J A > 1/2. The vertical and horizontal dot-dashed 
lines show the asymptotes for g —* J A and g/JA — > oo, re- 
spectively. The results refer to the range where the BP and 
LDP bands are far from each other, compared to the LDP 
bandwidth \J\. These bands are sketched in the inset (the 
BP band is above the LDP band for JA/g > 1). The dashed 
line in the inset shows where the energy level of the localized 
BP state is located with respect to the bands. 



(26) 



IV. ANTIRESONANT DECOUPLING OF 
TWO-EXCITATION STATES 



The analysis of the preceding section does not apply if 
the pair binding energy JA is close to the defect excess 
energy g. When \g — JA\ is of the order of the LDP 
bandwidth | J|, the BP and LDP states are in resonance, 
their bands overlap or nearly overlap with each other. 
One might expect that there would occur mixing of states 
of these two bands. In other words, a delocalized magnon 
in the LDP band might be scattered off the excitation on 
the defect, and as a result they both would move away as 
a bound pair. However, as we show, such mixing does not 
happen. In order to simplify notations we will assume in 
what follows that g,J,A > 0. 
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A. A bound pair localized next to the defect 

In the resonant region we should reconsider the analy- 
sis of the next-to-the-defect bound pair localized on sites 
(no + l,n + 2) or (n — 2,n — 1). As \g — JA| de- 
creases, the energy of this pair i|26[l moves away from the 
BP band, see Fig. [21 At the same time, the distance be- 
tween the pair energy and the LDP band E BP — 2e\ — g 
becomes smaller with decreasing \g — JA| as long as 
| JA - g\ > \ J\/2. For | JA - g\ ~ J, the next-to-the- 
defect pairs are hybridized with LDP's. The hybridiza- 
tion occurs in first order in the nearest-neighbor coupling 
constant J. 

To describe the hybridization, we will seek the solution 
of the Schrodinger equation (|10|) in the form of a linear 
superposition of an LDP state H18|) and a pair on the next 
to the defect sites, 

a(n,m) = (Ce i9m + C e' Wm )5 n ^ 

+ a(no + l,n + 2)6 nt n +iSm,n +2 

+ a(n Q - 2, n - l)<5„,„ -2<S m ,ri -i (27) 

(we remind that to > n). The energy of a state with 
given 9 can be found from (|10l) as before by considering 
to far away from no- It is given by Eq. 119|) . 

The interrelation between the coefficients C, C , et(no + 
1, no + 2), and a(no — 2, no — 1), as well as the values of 9 
should be obtained from the boundary conditions. These 
conditions follow from the fact that the energy of a pair 
on sites (no, no ±1) and the energies of the pairs de- 
scribed by Eq.(23 differ by - g » JA, cf. (gDJl. There- 
fore the pairs on sites (no, no ±1) are decoupled from 
the states l|27jl. and in the analysis of the LDP's we can 
set a(no,no ± 1) = 0. Decoupled also are unbound two- 
excitation states with no excitation on the defect, i.e., 
two-magnon states. Therefore a(n, m) = if simultane- 
ously m — n > 1 and (n — no) (to — no) =/= 0. 

From Eq. (|1U[) written for n = no,ni = no + 2 and 
n = no — 2, m = riQ [with account taken of the relation 
a(n, to) = a(m, n + N)] we obtain 

a(n + l,n + 2) = Ce ie{n ° +l) + c' e - i9 ( n °+i\ (28) 
a(n - 2, n - 1) = Ce ie( - no+N -^ + c > e -ie(n +N-i) _ 

The interrelation between C and C and the equation 
for 9 follow from Eqs. J2ZJ), and (JTUJ) written for 

n = no + 1, m = no +2 and n = no — 2, m = no — 1. They 
have the form 



a = 



2( JA — g)e ie — J 
2(JA- g)e~ ie - J 



Ce 



2i0n o 



and 



f(9) = f(-9), 
f(9) = e t9N [2(JA - g)e- w - j\ 



(29) 



(30) 



Equation (|3U[1 is an 2iVth order equation for exp(i6>). 
Its analysis is completely analogous to that of Eq. (|23|l . 



The roots 9 — 0, 7r are spurious, and the roots 9 of op- 
posite signs describe one and the same wave function. 
Therefore Eq. 1)30(1 has N — 1 physically distinct roots 9. 
Real 6*'s correspond to extended states. Complex roots 
appear for 2| JA — g\ > J. These roots, 0^ DP , describe 
localized states. In the limit of a long chain, N — > oo, 
the imaginary part of one of them is 



lm9 



W 

LDP 



]n\2(JA-g)/J\, 



(31) 



whereas the other root has just opposite sign. 

The wave function of the localized state is maximal 
either on sites (n + l,n + 2) or (n — 2,n — 1) and 
exponentially decays into the chain. For JA — g < this 
decay is accompanied by oscillations, Re 9^' DP = tt. The 
energy of the localized state is 



E 



LDP 



2ei + JA ■ 



(J/2) 2 
JA-g- 



(32) 



The localized state (|2Tjl . (|5ll is a "surface- type" state 
induced by the defect. It is the resonant-region analog 
of the localized next-to-the-defect state discussed in Sec- 
tion III.B3. The wave function of the latter state 122L 
(|25|l was a linear combination of the wave functions of 
bound pairs. In contrast, the state given by Eqs. i|2"7)l . 
(|31|l is a combination of the wave functions of the bound 
pair located next to the defect and a localized-delocalized 
pair. 

The evolution of the surface-type state is controlled 
by the difference between the excess energies of binding 
two excitations in a pair or localizing one of them on the 
defect | JA — g\. As | J A — g\ varies, the state changes in 
the following way. It first splits off the BP band when 
| JA — g | becomes less than g, see Fig. |21 Its energy 
moves away from the band of extended BP states with 
decreasing | JA — g\ and the localization length decreases 
[cf. Eq. lj2"51) ]. Well before | JA — g\ becomes of order J, 
the state becomes strongly localized on the sites (n + 
l,7io + 2) or (n - 2,n - 1). 

In the region | J A — g\ ~ J the localized state becomes 
stronger hybridized with LDP states than with extended 
BP states. This hybridization occurs in first order in J, 
via a transition (no + 1, no + 2) — > (no, no + 2) [or (no — 
2, no — 1) — ► (no — 2, no)]. In this region the localization 
length increases with decreasing | J A — g\, cf. Eq. 131fl . 
Ultimately, for |JA — g\ = J/2 the localized surface-type 
state disappears, as seen in Fig. |3J The evolution of the 
energy of the localized state with JA — g is shown in 
Fig. El 

The crossover from hybridization of the localized 
surface-type BP state with extended BP states to that 
with LDP states occurs in the region \g\ ^> | J A — g\ 3> 
| J|. It is described using a different approach in Ref.l20l 
We note that the expressions for the energy of the lo- 
calized states (|26|) and (I32|l go over into each other for 
|JA- 5 |«| 5 |. 



9 




2(JA-g)/J 

FIG. 3: The energy difference between the localized BP state 
and the LDP band, eldp = {E^bp - 2ei - g)/(J/2), vs the 
energy mismatch 2(JA — g)/J. In the region between the 
dot-dashed lines the BP state is delocalized. The results refer 
to the case where the BP and LDP bands are close to each 
other or overlap, as sketched in the inset (the BP band is 
above the LDP band for J A > g). The dashed line in the 
inset shows where the energy level of the localized BP state 
is located with respect to the bands. 

B. Decoupling of bound pairs and LDP's in the 
resonant region 

We are now in a position to consider resonant coupling 
between extended states of bound pairs and localized- 
delocalized pairs. To lowest order in A -1 , a BP on sites 
(n, n + 1) far away from the defect can resonantly hop 
only to the nearest pair of sites, see Eqs. (|15fl . 1)16(1. As 
noted before, the hopping requires an intermediate vir- 
tual transition of the BP into a dissociated state, which 
differs in energy by JA. 

A different situation occurs for the BP on the sites 
(ri + 2, n + 3) [or (n — 3,n — 2)]. Such BP can hop onto 
the sites (no + 1 , no + 2) and (no + 3, no + 4) , as described 
by (|16|) . But in addition, for | J A — g\ < J it can make 
a transition into the LDP state on sites (no, no + 3) [or 
(no — 3, no)] . Indeed, such state has the energy sa 2e\+g, 
which is close to the BP energy » 2ei + JA. 

The transition (no + 2, no + 3) — > (no, no + 3) goes 
through the intermediate dissociated state (n + l, no +3), 
which differs in energy by ~ JA. It can be taken into 
account by adding the term 5H^' to the BP hopping 
Hamiltonian Ijltifl for the sites (no + 2, no + 3), 

5H {t) a{no + 2, n + 3) = ( J/4A)a(n , n + 3). (33) 

Extended BP states are connected to LDP states only 
through a BP on the sites (no + 2, no + 3). We are now 
in a position to analyze this connection. From Eqs. Hlti|) 
and we see that 

[H® + 6H (t) }a(n + 2, n + 3) = ( J/4A) 

x [a(n , n + 3) + a(n + 1, n + 2)] + A, (34) 

where A is a linear combination of the amplitudes a (no + 
2, no + 3) and a(no + 3, no + 4). 



The sum a(n , n + 3) + a(n + l, n + 2) in Eq. ffi^i can 
be expressed in terms of the LDP wave functions 127fl . 
With account taken of the interrelation (|29|l between the 
coefficients C, C in Eq. (|27|l . we have 

a(n + 1, n + 2) + a(n , n + 3) 
ocCsin0cos0[JA-3- Jcos0]. (35) 

Equations l|34ll . 135fl describe the coupling between the 
BP on the sites (no + 2, no + 3) and the LDP cigenstates 
with given 9. 

An important conclusion can now be drawn regarding 
the behavior of BP and LDP states in the resonant re- 
gion. The center of the BP band lies at 2ei+ JA+( J/2A) 
(HU, and the BP band is parametrically narrower than 
the LDP band 2ei + g + J cos 9 (JTHJ . When the BP band 
is inside the LDP band, this means that, for an appropri- 
ate wave number 9 of the LDP magnon, JA = g + Jcos9, 
to zeroth order in A -1 (this is the approximation used to 
obtain the LDP dispersion law). It follows from Eqs. I|34|) 
and (|35|l that, for such JA — g and 9 there is no coupling 
between the LDP and extended BP states. 

The above result means that LDP and extended 
BP states do not experience resonant scattering into 
each other, even though it is allowed by the energy- 
conservation law. Such scattering would correspond to 
the scattering of a magnon off the excitation localized 
on the defect, with both of them becoming a bound pair 
that moves away from the defect, or an inverse process. 

Physically, the antiresonant decoupling of BP and LDP 
states is a result of strong mixing of a BP on the sites 
(no + 1, no + 2) and LDP's. Because of the mixing, the 
amplitudes of transitions of extended BP states to the 
sites (no + l,no + 2) and (no, no + 3) compensate each 
other, to lowest order in A -1 . 

To illustrate the antiresonant decoupling we show in 
Fig. 0] two types of time evolution of excitation pairs. In 
both cases the initial state of the system was chosen as a 
pair on the sites (no + 1, no + 2), i.e., a(no + 1, no + 2) = 1 
for t — 0. The time-dependent Schrodinger equation 
was then solved with the boundary condition that corre- 
sponds to a closed chain. 

The solid lines refer to the case of nonoverlapping BP 
and LDP bands, | J| <C \g - JA| and \g\ < | JA|/2. In 
this case there does not emerge a localized surface-type 
BP state next to the defect. Therefore an excitation 
pair placed initially on the sites (no + l,no + 2) reso- 
nantly transforms into extended BP states and propa- 
gates through the chain. This propagation is seen from 
the figure as oscillations of the return probability and the 
probability to find the BP on another arbitrarily chosen 
pair of neighboring sites (n + 2, n + 3). The oscillation 
rate should be small, of the order of the bandwidth J/ A. 
This estimate agrees with the numerical data. It is seen 
that the BP state is not transformed into LDP states. 
The amplitude of LDP states on sites (no,m > no + 2) 
remains extremely small, as illustrated for m = no + 2. 

The dotted lines in Fig. 0] show a completely different 
picture which arises when the BP band is inside the LDP 
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band. In this case an excitation pair placed initially on 
(no + 1, + 2) hybridizes with LDP rather than BP 
states. A transformation of the pair (no + 1, no + 2) into 
LDP's with increasing time is clearly seen. The period of 
oscillations is of the order of the reciprocal bandwidth of 
the LDP's J -1 , it is much shorter than in the previous 
case. Remarkably, as a consequence of the antiresonance, 
extended states of bound pairs are not excited to any 
appreciable extent, as seen from the amplitude of the 
pair on the sites (n + 2, n + 3). 
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FIG. 4: Time evolution of a two-excitation wave packet in 
an XXZ chain with a defect; \a(n, m)| 2 is the occupation of 
sites (n, m). Initially an excitation pair is placed next to the 
defect on sites (no + l, no + 2). The solid and dashed lines refer 
to the cases where the BP and LDP bands are, respectively, 
far away from each other (g — JA/4) and overlapping (g — 
J A) . In the first case a bound pair slowly oscillates between 
neighboring sites (n, n+1) and does not dissociate [a(no,no + 
2) remains very small]. In the case of overlapping bands, the 
pair at (no + l, no+2) is hybridized with LDP's, but practically 
does not mix with bound pairs on other sites [a(no + 2, no + 3) 
remains very small]. The results refer to a ten-site closed chain 
with A = 10. 



V. CONCLUSIONS 

We have analyzed the dynamics of a disordered spin 
chain with a strongly anisotropic coupling in a magnetic 
field. A defect in such a chain can lead to several localized 
states, depending on the number of excitations. This is 
a consequence of the interaction between excitations and 
its interplay with the disorder. We have studied chains 
with one and two excitations. 

The major results refer to the case of two excitations. 
Here, the physics is determined by the interrelation be- 
tween the excess on-site energy of the defect g and the 
anisotropic part of the exchange coupling J A. Strong 
anisotropy leads to binding of excitations into nearest- 
neighbor pairs that freely propagate in an ideal chain. 
Because of the defect, BP's can localize. A simple type of 
a localized BP is a pair with one of the excitations located 
on a defect. A less obvious localized state corresponds to 
a pair localized next to a defect. It reminds of a sur- 
face state split off from the band of extended BP states, 
with the surface being the defect site. We specified the 
conditions where the localization occurs and found the 
characteristics of the localized states. 

Our most unexpected observation is the antireso- 
nant decoupling of extended BP states from localized- 
delocalized pairs. The LDP's arc formed by one excita- 
tion on the defect site and another in an extended state. 
The antiresonance occurs for g w JA, when the BP and 
LDP bands overlap. It results from destructive quan- 
tum interference of the amplitudes of transitions of BP's 
into two types of resonant two-excitation states: one is an 
LDP, and the other is an excitation pair on the sites next 
to the defect. As a result of the antiresonant decoupling, 
extended BP's and LDP's do not scatter into each other, 
even though the scattering is allowed by energy conser- 
vation. This means that an excitation localized on the 
defect does not delocalize as a result of coupling to other 
excitations. 

The occurrence of multiple localized states in the pres- 
ence of other excitations is important for quantum com- 
puting. It shows that, even where the interaction be- 
tween the qubits is "on" all the time, we may still have 
well-defined states of individual qubits that can be ad- 
dressed and controlled. One can prepare entangled local- 
ized pairs of excitations, as we discussed in Sec. Ill B 2, or 
more complicated entangled excitation complexes. The 
results of the paper also provide an example of new many- 
body effects that can be studied using quantum comput- 
ers with individually controlled qubit transition energies. 



Acknowledgments 

This research was supported in part by the NSF 
through Grant No. ITR-0085922 and by the Institute 
for Quantum Sciences at Michigan State University. 



11 



APPENDIX A: ONE EXCITATION 

In an infinite spin chain in a magnetic field, the 
anisotropy of the spin-spin interaction does not affect the 
spectrum and wave functions of one excitation. The ma- 
trix element of the term ^ n cr^n+i m the Hamiltonian 
||TJ is just a constant. However, the situation becomes 
different for a chain of finite length, because the cou- 
pling anisotropy can lead to surface states. In the case 
of two excitations, analogs of surface-type states emerge 
near defects in an infinite chain, as discussed in Sec. III. 
Here, for completeness and keeping in mind a reader with 
the background in quantum computing, we briefly out- 
line the results of the standard analysis of a finite-length 
spin chain with one excitation. 



1. Localized surface and defect-induced states in 
an open chain 

The Schrodinger equation for an excitation in an open 
chain has the form JHJ, and its solution a(n) to the left 
and to the right from the defect can be written in the form 
of a superposition of counterpropagating plane waves, 
Eq. ijfjj. The relation between the amplitudes of these 
waves C[ and C; lT - follows from the boundary conditions 
o(0) = a(N+l) = 0. By substituting Eq. © into Eq. © 
with n = 1 and n = N, we obtain 



D -- 



- -DQ, C' r = -D- l e 2i ^ N+ ^C r , (Al) 
[1 - Aexp(«0)]/[1 - Aexp(-i0)]. 



The relations between C! r , Ci jT , and the amplitude of 
the wave function on the defect site a (no) follow from 
Eqs. |J5| and JBJ for n = no ± 1, 

a(n ) = Cie i9no [l - De- 2l9n °] 

= C r e i9no [l - £)-l e 2^(iV+l-no)l _ ( A2 ) 

With I|A1(I and (|A2(I , all coefficients in the wave function 
© are expressed in terms of one number, a(no). It can 
be obtained from normalization. 

In a finite chain, the values of 9 are quantized. They 
can be found from Eq. JSJ with n = n . With account 
taken of (| A2|) . this equation can be written as 



2no + l)],(A3) 



f N {9) - D 2 f N (-6) = -{ig/JsmO) 
x lf N (6) + D 2 f N (-8) - 2D cos 9(N 



where 



f N (6)=exp\i6(N + l)]. 



(A4) 



The analysis of the roots of Eq. I|A3|) is standard. This 
is a 2(iV+ l)-order equation for exp(i0), but its solutions 
for 9 come in pairs 9 and —9. Each pair gives one wave 
function, as seen from Eq. ©. In addition, Eq. I|A3(I has 
roots 9 = 0, 7r; they are spurious (unless | A| = 1 and g = 



0) and appear as a result of algebraic transformations. 
Therefore Eq. I|A3[) has N physically distinct roots, as 
expected for a chain of TV spins. We note that, for g = 
the position of the impurity no drops out from Eq. I|A3() . 
and then the equation goes over into the result for an 
ideal chain. 

Solutions of Eq. (|A3() with real 9 correspond, in the 
case of a long chain, to delocalized magnon-type excita- 
tions propagating in the chain. Their bandwidth is 2\J\. 

Along with delocalized states, Eq. i|A3(l describes also 
localized states with complex 9. Complex roots of 
Eq. (|A3|) can be found for a long chain, where |Im#|(iV — 
n ), |Im6*|n » 1. They describe surface states, which 
are localized on the chain boundaries, and a state local- 
ized on the defect. The localization length of the states 
is given by |l/lm0|. 

The surface states arise only for the anisotropy param- 
eter |A| > 1. The corresponding values of 9 are 



9 S = ±iln|A| +7r6(-A). 



(A5) 



Here, the signs + and — refer to the states localized on 
the left and right boundaries, respectively, and Q(x) is 
the step function. 

From Eq. @ , the energy of the surface state is 



E s =£1 + J(A 2 + 1)/2A 



(A6) 



It lies outside the energy band of delocalized excitations. 
We note that, for A > 1 the surface states decay mono- 
tonically with the distance from the boundary (Re 9 = 0). 
For sufficiently large negative A, on the other hand, the 
decay of the wave function is accompanied by oscillations, 
and a(n) changes sign from site to site. 

A defect in a long chain gives rise to a localized one- 
spin excitation for an arbitrary excess energy g [Isj . The 
amplitude a(n) decays away from the defect as 

a(n) = a(n Q ) exp(i9 d \n - n |), (A7) 
9 d = tsmh- 1 (\g/J\) + Tre(-g/J). 

The energy of the localized state is 

E d = e 1 + (g 2 + J 2 ) 1/2 sgn 5 . (A8) 

For small \g/J\, we have Im 9 d ~ \g/J\, he., the re- 
ciprocal localization length is simply proportional to the 
defect excess energy \g\. In the opposite case of large 
\g/J\ we have Im 9 d = In |2<jr/J|. In this case the am- 
plitude of the localized state rapidly falls off with the 
distance from the defect, a(n) cx (2\g/ J|)~l"~"ol . 

If the localization length is comparable to the chain 
length, the notion of localization is not well defined. 
However, when discussing numerical results, one can for- 
mally call a state localized if its wave function exponen- 
tially decays away from the defect and is described by a 
solution of Eq. (|A3|) with complex 9. This is equivalent 
to the statement that the state energy E d lies outside the 
band of magnons in the infinite chain. In an open finite 
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chain such localized state may emerge provided the local- 
ization length is smaller than the distance from the defect 
to the boundaries. This means that the defect excess en- 
ergy \g\ should exceed a minimal value that depends on 
the size of the chain. The comparison of Eq. (|A7|) with 
the numerical solutions of the full equation I|A3J| for a 
finite chain is shown in Fig. [5J 

1.5 i ■ , ■ , ■ 1 




FIG. 5: The reciprocal localization length Im Od for an infinite 
chain [Eq. l|A7(l ] as a function of the defect excess energy g 
for J — 1 (solid line). Also shown are the results for an open 
chain with the number of sites N — 6 and 12 and no = N/2 
(stars and crosses, respectively) and a closed chain with 6 and 
12 sites (circles and squares, respectively). They are obtained 
from Eq. igU with A = 10 and from Eq. (jATTjl . In an open 
chain, solutions with nonzero Im 6d emerge starting with a 
certain \g/J\ > \g/J\mm- In a closed chain with even N there 
is no threshold in \g/J\ for the onset of states with Im 6d ^ 0. 



2. One-excitation states in a closed chain 

As pointed out in Sec. II, for a closed chain the solu- 
tion of the Schrodinger equation can be also sought in the 
form of counterpropagating waves with different ampli- 
tudes, Eq. Q. Clearly, the phases 9 and —9 describe one 
and the same wave function. The one-excitation energy 
E\ is given by Eq. JSJ. 

The interrelation between the amplitudes of the waves 
C, C in Eq. (7J and the amplitude of the wave function 
on the defect site a (no) can be obtained from Eq. (JSJ) 
with n = tiq ± 1. This equation has two solutions, 

e l9N = 1, a(n ) = Ce i9n ° + C'e~ l9no (A9) 

and 

a{n ) = Ce l0n ° (l + e l0N ) (A10) 



= C'e~' l9no (1 + e' WN ) [exp(i6N) ^ 1]. 
In order to fully determine the wave function a(n) 

0, Eqs. (|A"9)) . IjAlOjl should be substituted into the 
Schrodinger equation (JSJ for n = uq. 

Equations @ and (| A9|) can be satisfied provided that 
either g = 0, which means that there is no defect, or 
a (no) = 0. The first condition describes excitations in an 
ideal closed chain and is not interesting for the present 
paper. The condition a(rio) = corresponds to the wave 
function a(n) oc sin#(n — no), which has a simple physical 
meaning. It is a standing wave in an ideal chain with a 
node at the location of a defect. Because of the node, the 
corresponding state "does not know" about the defect, 
and therefore it is exactly the same as in an ideal chain. 

In the presence of a defect, the solutions of Eq. (|A9|> 
in the range of interest < 9 < ir are 9 = 2irk /N with 
k = l,2,...,(N-l)/2for oddN, orfc = 1, 2, . . . , JV/2-1 
for even N. 

The equation for 9 that follows from Eqs. © and (|A1 0|) 
has the form 

exp(i9N) -1 = - 19 [exp(i6>AT) + 1]. (All) 
J sin 9 

For g this equation has either (N + l)/2 (for odd N) 
or N/2 + 1 (for even N) solutions for ±9. Therefore the 
total number of solutions for 9 that follow from Eqs. (|A9|) 
and IjAllfl is AT, as expected. 

By rewriting Eq. (fATTJ) as tan(6W/2) = -(g/Jsin9) 
and plotting the left- and right-hand sides as functions 
of 9 (cf. Ref.[l]|, one can see that all physically distinct 
roots of this equation but one are real and lie in the 
interval < 9 < tt [except for one case, see below]. Such 
solutions describe delocalized states with sinusoidal wave 
functions. 

The complex root of Eq. (|A11|1 . 9 = 64, describes a 
state localized on the defect. For a long chain, Im 9^N ^> 

1, the solution has the form l|A7(l . as expected. An in- 
teresting situation occurs for a shorter chain. If N is 
even or if g/ J > 0, a localized solution with complex 9d 
emerges for any defect excess energy g. Thresholdless lo- 
calization does not happen in an open chain. In a closed 
chain, it arises because there is no reflection from bound- 
aries. For small positive g/ J one obtains the complex 
solution of Eq. l|AlT|l in the form 9 d w i(2g /N J) 1 /' 2 . The 
square-root dependence of Im 9d on g is seen from Fig.0 

Equation (jAllfl has a complex solution also for 
even N and small negative gj J. In this case 9d w 
i(—2g/NJ) 1 / 2 + 7T, i.e., the decay of the wave function 
a(n) is accompanied by sign flips, a(n + l)/a(n) < 0. 
Such oscillations cannot be reconciled with the period- 
icity condition for odd N. Therefore, for odd N and 
negative gj J a decaying solution arises only when —gj J 
exceeds a threshold value. One can show from Eq. (|A11|) 
that this value is \g/J\ m in — 2 /AT, which has also been 
confirmed numerically. 
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